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We present simulation results of the (2 + l)d four-fermion model with a baryon 
chemical potential ^. We examine temporal correlation functions of the vector me- 
son, and find evidence of phonon-like behavior characterised by a linear dispersion 
relation in the long wavelength limit. We also discuss the consistency of our nu- 
merical results with analytical solutions to the Boltzmann equation corresponding 
to zero sound. We argue that our results provide the first evidence for a collective 
excitation in a lattice simulation. 

We study the (2 -I- l)d Gross-Neveu (GN) model with a chiral sym- 
metry describing Nf flavors of self-interacting fermions, with continuum 
lagrangian 

= +^^o)V-~^(^V)2. (1) 

As /X is increased, the model exhibits a sharp first-order transition at a 
critical from a phase where a fermion mass is dynamically generated 
but baryon density ub = ('07o'0) vanishes to one where rrif = and ns (x 

[1]. The main new feature emerging for > /ic is the existence of a 
new scale, the Fermi momentum kp. Since the lowest energy excitations 
of the ground state have « kp, measurements of Euclidean timeslice 
correlators with k ^ are mandatory. 

If certain conditions are fulfilled, such as an effective interaction be- 
tween quasiparticles which is both short-ranged and repulsive, then there 
is a massless bosonic excitation in the spectrum as T — > 0. Carrying zero 
baryon number, it is a phonon, ie. a quantum of a collective excitation 
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called zero sound [2]. Sound propagation occurs in any clastic medium; 
zero sound happens when the elasticity originates not from collisions be- 
tween individual particles, but from the force on a single particle due to its 
coherent interaction with all others present in the medium. It can be pic- 
tured as a propagating distortion in the local shape of the Fermi surface, 
and its speed of propagation exceeds that of conventional "first" sound. 
In [3] we solved the self-consistent equation for zero sound in the frame- 
work of Fermi liquid theory in the case of the id GN model to leading 
non-trivial order in 1/Nf. We found solutions for sound propagation with 
speed /3o > 0f, the Fermi velocity. Here we discuss our numerical results 
and compare them with the analytical solutions presented in detail in [3]. 

The lattice regularized action of our model uses staggered fermions. 
Details of the formulation and simulation algorithm are given in [1,4]. The 
simulation parameters are Nf = A and a/g'^ = 0.75, corresponding to a 
physical fermion mass at /i = of m/ ~ 0.17a"^. All results are taken in 
the chirally-restored phase, ie. with > jjc — 0.16a~^. 

By analysing the decay of the appropriate correlation function with 
Euclidean time, we calculated the dispersion relations E{k) for the spin-^ 
quasiparticle, which carries a baryon charge, and co{k) for various meson 
states of the form 'tpTtp which are best thought of as excitations formed 
from a particle-hole pair. We simulated x 48 systems with Lg = 32 
and 48, with fi ranging from 0.2 to 0.6, and measured E{k),Lu{k) for k = 
(fc,0) with k = 0, 27r/Ls, . . . , 7r/2. The increased momentum resolution 
offered by Lg = 48 has proved to be important. Note that the staggered 
fermion action is only invariant under translations of two lattice spacings, 
which restricts the space of accessible momenta. In [3,5] we measmed the 
quasiparticle propagator and showed that the transition between particle 
and hole behavior is rather sudden, characteristic of a well-defined Fermi 
surface. We also presented [5] the theoretical predictions of Fermi liquid 
properties to leading non-trivial order in 1/iV/ for /z ^ /Uc. The results for 
the Fermi velocity Pp and the Fermi momentum kp are given by: 

As discussed in [3,5] the agreement between the analytical predictions and 
the numerical results on 48^ lattices for /i = 0.2, 0.6 and with discretiza- 
tion effects taken into account is at best qiialitativc. It should be noted that 
analytic predictions of Pf and kp to 0{1/Nf) for finite /j. strictly requires 
a knowledge of the gap equation to two loops, which is not yet available. 
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For the remaining analysis we will therefore assume the free field values 
kp = ix, I3f = 1- 

Next we consider the meson sector, ie. correlators of the form 

Cr(fc, t) = 5^(v;rv(o, 0) virv(f , t))e^'^-^ (3) 

X 

which carry zero baryon number. In [5] we showed that the generic 
large-distance behaviour in any given channel is dominated by zero-energy 
particle- hole pairs, and as a result the decay is algebraic, ie. Cr(t) oc t~^^^\ 
with the exponent A determined by the geometry of the overlap between two 
Fermi disks with relative displacement k between their centres. A particu- 
larly favourable configuration occurs for |fc| » 2/Lt, in which case the Fermi 
surfaces just kiss, and A = |. In this paper wc focus on point-split me- 
son operators corresponding to spatial components of the conserved vector 
current, ie. of the form 

J^ {x) = ^ [x{x)x{x + i) + x{x + i)x{x)] , (4) 

where XtX ^re staggered fermion fields, rji{x) = (—1)*, and rj2{x) = 
(— 1)*+"^!. We set k = (fc,0), and define Cy in terms of the correlator 
0'i(0)ii(a;)) and Ci_ ~ (i2(0)i2(a;)). In Fig. 1 we show data 
taken on a 48'^ lattice at /i = 0.5, for both k — ^ 2^ and k = tt — ^ 
(the modulus is taken because we use a log-linear scale, and C fluctuates 
in sign). Whereas C± and C\\{k = |) all show behaviour consistent with 
algebraic decay, the decay in the C||(fc = n — ^) channel is much faster, 
and resembles the exponential decay expected of an isolated pole. We have 
fitted it with the form 

q|(7r-fc,t) = Aexp{~-n{k)t) + {-IYB cxp{-u;{k)t), (5) 

in most cases employing datapoints with t G [10, 38]. For small k, the coef- 
ficient B ^ A, suggesting that the correlator is dominated by a pole in the 
alternating channel. In Fig. 1 wc plot the resulting dispersion relation 
for /z = 0.2, . . . , 0.6. The behaviour w oc A: as A: ^ suggests the presence of 
a massless pole similar to that of a phonon. Support for this interpretation 
comes from recasting the meson bilinear in terms of continuum-like fields 
q"°', q""" having spinor index a = 1, . . . , 4 and 'color' index a = 1, 2 [6]. We 
obtain 

{-in-lYx{x)VlX{x + l)-tq{-fo<SiT;)q, (6) 

demonstrating that the excitation can be viewed as an oscillation of local 
baryon density. 
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Figure 1. Left: Meson propagators |C|| | and \C±_ \ for both k = ^ and ti" — for ^ = 0.5. 
The dashed line shows a fit of the form (5). Right: Dispersion relation u}{k) for various 
11 for the channel defined by the correlator Cy (tt — fc). Data are from 48^. 



Although we are examining the k ^ Q hmit, we are trying to model a 
Fermi surface phenomenon on the lattice; the appropriate scale with which 
to compare (3^ is the "bare" Fermi velocity (ie. with no discretisation 
correction) which for free fields with E{k) = — /i + sinh^^ (sin fc) is given by 



Pf — 



dE 
dk 



1-sinhV^ (7) 



E=0 \ ^ + M 

The general trend in our data is that the speed ratio s — / (i^'"^ increases 
towards unity as n increases (for ji = 0.2 s = 0.896(8) and for fi = 0.6 
s ~ 0.959(2)). In [3] we used Fermi liquid interaction to leading non- 
trivial order in 1/Nf to find analytical solutions to the Boltzmann equation 
corresponding to zero sound. The excitation is spin and isospin symmetric, 
and has speed ratio s > 1 for almost all fi > fic- To what extent can we be 
sure that our numerical results describe the same physical phenomenon? 
Solutions with s < 1 allow the possibility of emission of a phonon with 
momentum q and energy (3o\q\ from a quasiparticle with momentum k and 
energy fi + (3p{\k\ — kp). This is allowed kinematically for s < 1, the angle 
of emission cj) satisfying 

cos(f) = s + ^(1- s^) > s. (8) 
2|fc| 

All radiation is emitted within a cone of half- angle cos^^ s centred on the 
quasiparticle trajectory. This well-known phenomenon is variously known 
as Landau damping, Cerenkov radiation, or most appropriately in the cur- 
rent context, as a sonic boom. 
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First let us argue how our numerical results can support a state resem- 
bling a simple pole but with s < 1, in apparent contradiction to the above. 
On a spacetime lattice, the radiation process is constrained because there is 
a natural lower boimd for the angle of emission, ^min ~ 27r/Ls|9|- Landau 
damping is thus kinematically forbidden for 

S > COS(^iniii ^ 1 - roi^o - (9) 

With a conservatively high |g| ~ /i/2, then on a 48^ lattice wo have no 
damping for s > 0.73 at ^ = 0.2, rising to s > 0.97 at /i = 0.6. It is 
plausible therefore that Landau damping is suppressed in our lattice data 
and that the phonon is described by an isolated pole even if s < 1. 

Simulations on volumes considerably larger than those used here, how- 
ever, will be needed to disentangle the various possible systematic effects 
due to finite Ls, finite n/T, and non-zero lattice spacing in order to deter- 
mine the sign of s — 1 for finite Nf. We think the most probable explanation 
for ,s < 1 is an effect of non-zero T = (aLt)^^; the physical temperature 
of the Lt = 48 lattice decreases by a factor of roughly two going from 
g~'^ = 0.75 to g~'^ = 0.6, and according to our lattice data |s — 1| decreases 
by roughly the same factor. It is known empirically that the speed of sound 
in liquid ^He increases as T — » and the dominant mode of propagation 
changes from first sound to zero sound [7]. In principle, since our extracted 
value of s i$ 1 has relied on a rescaling of /3f to take account of discreti- 
sation artifacts, we also need to go to considerably finer lattices, ie with 
^baro ~ X, in order to demonstrate a clear distinction between our signal 
and first sound with expected propagation speed (3i ~ 0.7. If, however, 
the theoretical arguments about zero sound dominating as T ^ can be 
taken seriously, then for the first time we have succeeded in identifying a 
collective oscillation in a lattice simulation. 
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